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We compute non-perturbatively the renormalization constants of composite opera- 
tors for overlap fcrmions by using the regularization independent scheme. The scaling 
behavior of the renormalization constants is investigated using the data from three lat- 
tices with similar physical volumes and different lattice spacings. The approach of the 
renormalization constants to the continuum limit is explored. 
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1. Introduction 

Following previous papers ^2 i n which we computed non-perturbatively the renor- 
malization constants of composite operators with overlap fermions in quenched 
QCD, this paper will study the scaling behavior of the renormalization constants. 

We adopt the non-perturbative renormalization method which was introduced 
by Martinelli et al. ^ The method allows a full non-perturbative computation of 
the matrix elements of composite operators in the Regularization Independent (RI) 
scheme^S (it is called the RI' scheme by ChetyrkinEu). The matching between the 
RI scheme and MS, which is intrinsically perturbative, is computed using continuum 
perturbation theory, whi ch i s well behaved. 

The overlap fermion was proposed by Narayanan and Neuberger to evade 
the so called "no-go" theorem^. The action in the massless limit pres erve s a lattice 
form of chiral symmetry even at finite lattice spacing and volume lil^. The use 
of the overlap action entails many theoretical advantages: it has no additiv e m ass 
renormalization, there are no order a artifacts, and it has very good scaling 

2. Non-perturbative renormalization method 

The renormalized operator 0(/i) is related to the bare operator, 0(a), calculated 
on the lattice via 

0(jx) =Zo(m,9(a)) 0(a), (1) 
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Table 1. Lattice parameters. 



Action 


Size 


Nsamp 


P 


a (fm) 


uo 


Physical Volume (fm 4 ) 


Improved 


16 3 X 32 


500 


4.80 


0.093 


0.89650 


1.5 3 X 3.0 


Improved 


12 3 X 24 


500 


4.60 


0.123 


0.88888 


1.5 3 x 3.0 


Improved 


8 3 x 16 


500 


4.286 


0.190 


0.87209 


1.5 3 x 3.0 



In this work, we will consider the fermion operators 

O r (x)=^(x)T^(x), (2) 

where T are the Dirac gamma matrices T £ {1, 7 M , 75, 7^75, u^} and the corre- 
sponding notations will be { S, V, P, A, T } respectively. 

The renormalization condition is imposed directly on the three-point vertex 
function To{pa), which is calculated in a fixed gauge, e.g., Landau gauge in our 
case, at a momentum scale p 2 = p 2 

r t \ 1 Z (pa,g(a)) 

tO,ren{Pa)\p2 =fl 2 = — — — -L {pa) p 2 = 2 = 1. (3) 

Z^itia, g(a)) 

1/2 

Here Z^ is the field or wave-function renormalization constant, ^ren = Z^ $ . 

It is apparent that we can only get the ratio of the renormalization constant 
Zq for the operator O and the wave-function renormalization constant Z^, from 
the renormalization condition of Eq. J^J. In order to obtain the renormalization 
constant Zq for the operator O, one needs to know Z^ first. In this work, we will 
obtain Z^ directly from the quark propagator. It can be defined from the Ward 



:tly 

M 



Identity (WI) as I 

1 Tr Eu=1.4 7/x(PM a ) 5 '(P a )" 



Z' 



-1 



12 4£_ M (p, 



a 



(4) 



which, in Landau gauge, differs from Z^ by a finite term of order a 2 . The matching 
coefficients have been computed using continuum perturbation theory El. 



3. Numerical details 

We work on three lattices, each with a different lattice spacing, a, but having simi- 
lar physical volume. Lattice parameters are summarized in Table 1. The quenched 
gauge configurations are created using a tadpole improved plaquette plus rectangle 
(Liischer-Weisz) gauge action and gauge fixed to the Landau gauge using a Conju- 
gate Gradient Fourier Acceleration algorithm. 
The overlap-Dirac operator we use is 

D( mq )=p+^ + (p- 1 ^h 5 e(H). (5) 

Where p is regulator mass and m q is bare quark mass, and e(H) is the matrix sign 
function of an Hcrmitian operator H — j^Dw- We use the tadple improved Wilson 
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kernel Dw in the overlap operator, and k — 0.19163 is used for the regulator mass 
p for all three lattices. We calculate the overlap quark propagator for 15 bare quark 
masses on three lattices, they are 53, 59, 71, 83, 94,106, 124, 142, 177, 212, 266, 
554, 442, 531, and 620 MeV respectively 

Our calculation begins with the evaluation of the inverse of the overlap-Dirac 
operator. After we calculate the quark propagator in coordinate space for each 
configuration, we use the Landau gauge fixing transformation matrix to rotate the 
quark propagator to Landau gauge. Then the discrete Fourier transformation is used 
to obtain the quark propagator in momentum space. Afterward, we calculate five 
projected vertex functions To{pa) for each bare quark mass, and extrapolate to the 
chiral limit. These projected vertex functions To(pa) are in general dependent on 
(pa) 2 . The dependence may come from two sources. One is from the usual running 
of the renormalization constant in the RI scheme. The other is from possible (pa) 2 
errors and we need to remove this. In order to confront experiment, it is preferable 
to quote the final results in the MS scheme at a certain scale. One needs to transform 
the resu lts i n the RI scheme to the MS scheme. The detailed analysis can be found 



4. Scaling behaviors 

We work on three lattices with similar physical volumes and different spacings a 
to investigate the the scaling behavior of the renormalization constants. Here we 
compare the results of renormalization constants Z^, Zy, Z$ and Zt on the dif- 
ferent lattices in the MS scheme at 2.0 GeV. Fig. 01 shows the four renormalization 
constants Z^, Zy, Zg and Zj< against the square of the lattice spacing a. Because 
overlap fermions are free of 0(a 2 ) errors, the leading term must be proportional to 
a 2 . We use a simple linear fit Zq = c\ + c 2 a 2 , and take C\ as the value of the renor- 
malization constant Zo in the continuum limit. The numerical values of calculated 
renormalization constants in the continuum limit in the MS and RI schemes at 2.0 
GeV are displayed in Table 2. 

In the continuum, the vector current is conserved, so Zy should be equal to one. 
For the axial vector current, due to the PCAC relation, it will be conserved at large 
momenta, where the anomaly has no effect. Our result for Zy and Za in Table 2 
compares favorably with 1. 



in Refs. GE1. 



Table 2. Results for Z in the continuum limit. 



Z - factor RI scheme at 2 GeV MS scheme at 2 GeV 




1.059±0.008 
0.987±0.011 
0.765±0.021 
1.069±0.010 



1.045±0.008 
0.987±0.011 
0.898±0.025 
1.047±0.010 
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Fig. 1. The renormalization constants Z^,, Z\, Z$, and in MS at 2.0 GeV against the square 
of the lattice spacing a. The straight line is the linear fit Zo = ci + C2a 2 . 
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